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ABSTRACT 

Using the NICMOS instrument on the Hubble Space Telescope, we have mea- 
sured the dayside spectrum of HD 209458b between 1.5-2.5 /xm. The emergent 
spectrum is dominated by features due to the presence of methane (CH 4 ) and wa- 
ter vapor (H2O), with smaller contributions from carbon dioxide (CO2). Combin- 
ing this near-infrared spectrum with existing mid-infrared measurements shows 
the existence of a temperature inversion and confirms the interpretation of pre- 
vious photometry measurements. We find a family of plausible solutions for the 
molecular abundance and detailed temperature profile. Observationally resolving 
the ambiguity between abundance and temperature requires either (1) improved 
wavelength coverage or spectral resolution of the dayside emission spectrum, or 
(2) a transmission spectrum where abundance determinations are less sensitive 
to the temperature structure. 
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Subject headings: planetary systems — techniques: spectroscopy 



Introduction 



HP 209458b has the d istinction of being the first identified transiting extrasolar planet 

( jCharbonneau et al. II2000I ). This 0.69 Jupiter-mass planet orbits a GOV stella r primary 

(K m „ g =6.3) every 3.52 days and was the first exoplanet with a detected atmosphere (jCharbonneau et al. 



20021 ). I t was also one of the first planets found to emit detectable amounts of infrared ra- 
diation (jDeming et al. II2005I ). The dayside mid-infrar ed spectrum for this planet, based on 



Spitze r IRS observations, has been previousl y obtained (IRichardson et al. ll2007l ; ISwain et al. 



2008a ). Mid- infrared dayside photo metry (Knutson et al. 20081 ) provides evidence for an 
atmospheric temperature inversion (IBurrows. et al. 1120071) and has been use d to support a 
classification scheme for the hot- Jupiter type exoplanets (IFortney et al. 1120081) . Mid-infrared 



phot ometric measurements of the orbital phase light curve of HD 209458b (j Cowan, et al. 



20071 ) find that the dayside and nightside brightness temperatures are similar. Models of 
atmosph eric circulation for HP 209458b suggest that the zonal wind structure i s a function of 
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altitude (IShowman et al. 1120081). a nd mid-infrared dayside emission could vary (IRauscher et al. 



20081 ). 



The recent announcement of the detection of H 2 and CH 4 in the terminator regions 
of the hot- Jupiter HP 189733b (Swain, Vasisht & Tinetti 2008b; hereafter SVT08) via near- 
infrared spectroscopy using the Hubble Space Telescope (HST) demonstrated that molecular 
spectroscopy of bright transiting planets is possible. The subsequent HST detection of 
H 2 0, C0 2 , and carbon monoxide (CO) on the dayside of HP 189733b provides significant 
constraints for atmospheric models and suggests that exoplanet molecular spectroscopy may 
become routine (Swain et al. 2009; hereafter S09) in the near future. This is significant 
because molecules provide a powerful tool for determining the conditions, composition, and 
chemistry of exoplanet atmospheres. In this paper, we report on the first ever spectroscopic 
detection of molecules in the atmosphere of HP 209458b, making this the second exoplanet 
for which molecular spectroscopy has now been demonstrated. 



2. Observations & Data Analysis 

We observed HP 209458 for five consecutive HST orbits using the NICMOS camera 
in imaging-spectroscopy mode covering the wavelengths 1.5-2.5 /im using the G206 grism. 
As in previous observations (SVT08, S09), the PEFOCUS mode was used, resulting in a 
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spectral resolution of R ~ 40. Observations were conducted on June 15, 2008, from UT 
09:38:59 to UT 16:45:24. A total of 310 usable snapshot spectra were acquired with an 



individual exposure time of T = 4.06 s. Using the system ephemeris (IKnutson et al. 1 120071 ) 
we determine that the bulk of the first two orbits (0\ & O2) coincides with the pre-ingress 
light curve, the third and part of the fourth orbit (O3 & O4) cover the occultation, while the 
fifth orbit (O5) is post-egress (see Fig. 1). We find that 75 exposures cover the full eclipse, 
43 exposures are in ingress/egress, and 192 exposures cover the pre-ingress or post-egress 
intervals. The effective exposure time for each spectrum was T = 4.06 s. To determine the 
spectrometer wavelength calibration, narrow band (F190) filter calibration exposures were 
acquired during 0\. 

Our data analysis methods, based on decorrelation of instrument parameters from the 
measured spectrophotometric light curves, have been described previously (SVT08, S09). 
For the results presented here, the data-analysis approach remains largely unchanged but is 
updated in two significant ways. First, the decorrelation for instrument parameters is based 
on detector pixel values; in the past we averaged the detector data to the defocused spectral 
resolution prior to the decorrelation. The new approach has the advantage of determining 
the corrected, decorrelated spectrum using the full oversampling of the instrument. Second, 
we now derotate and shift the spectrum (via interpolation) to achieve the same sampling at 
all wavelengths and times. This effectively places observations in a common, pixel-based, 
reference frame. We tested this new reduction method by computing a dayside spectrum of 
HD 189733b and comparing it to our previous result (S09); the agreement between the two 
methods was found to be excellent, with differences between the spectra always within the 
statistical uncertainty of the measurement. As with our previous observations of the primary 
and secondary eclipse of HD 189733b, we found the instrument parameters for the first orbit, 
0\, differed significantly from the remaining orbits, O2-O5; following our previous approach 
(SVT08, S09), we omitted the Oi spectra from the analysis. To add an additional constraint 
to any possible orbit-to-orbit brightness changes, we arranged for contemporaneous, high- 
precision j_vjsil2k_r2lwtom^try using the Micro-variability and Oscillations of Stars (MOST) 



satellite (IWalker et al. II2003I ). The MOST observations were timed to provide a photometry 



measurement during each HST orbit (se e Fig. 1), and the data were calibrated using the 



method described by iRowe et al. I (120081 ) 



3. Discussion 



The near-infrared dayside spectrum of HD 209458b shows three clear peaks at 1.6, 1.8- 
1.9, and 2.2 /j,m, while between the peaks the flux density decreases to near-zero (within the 
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measurement error). The spectrum has an average uncertainty per data point of ~ 7.5 x 10~ 5 , 
corresponding to a dynamic range of 13,300:1 (75 ppm), with a wavelength resolution of 
= 35 at 2 /mi. To maximize wavelength cove rage for the spectral interpretation, we 
incorporated mid-infrared Spitzer photometry data (IKnutson et al. 1120081 ) and spectroscopy 
data (originally obtained by Richardson e t al. 2007). We hav e previously undertaken our 
own calibration of these Spitzer IRS data (ISwain et al. Il2008a| ). and we use those results in 
the discussion that follows. 

The interpretation of the data was done using two independent spectral r e trieval meth- 
ods inc orporating; a line- by-line radiative transfer model ( [Griffith, et al. I (119981 ) : iTinetti et al. 
( l2007al lbL SVT08, S09), and the results from both models are in agreement. The radiative 
transfer calculations assume local thermal equilibrium (LTE) conditions, as expected for 
pressures exceeding 10~ 3 bar that are probed by the near-infrared spectrum. We evaluated 
a variety of temperature (T) together with the effects of CH 4 , H 2 a nd CO?, which are 



assum ed to have constant mixing ratios. We use the CH 4 line list of iNassar &: Bernath 



(120031 ) for temperatures at 800, 1000, and 1273 K to quantify the affects o f CH 4 at wave- 
l ength s of 1.58-2.5 /mi, outside o f which we use abso rption parameters by iRotham et al. . 
( 2005 ). The hot water lines from Barber et al. ( 2006 ) and Zobov et al. ( 20081 ) were used 



to quant ify the water features . Absorption by C0 2 is calculated from the HITRAN hot C0 2 
line list ( ITashkun et al. II2003I ). Absorption coefficients of H 2 and C0 2 are calculated using 
line-by-line techniques every 0.004 cm -1 wavenumbers. The radiative transfer calculation 
has 80 vertical grid points that extend from 10 -7 to 10 bar. The effects of particulates are 
excluded. 

The combination of near-infrared and mid-infrared data probes the temperature profile 
over a pressure range of a few bars to 10~ 5 bars and provides evidence for a temperature 
inversion. The 8 /im brightness temperature is 1770 K, while the 2.2-2.4 /im brightness 
temperature ranges from 1100-1600 K; this temperature difference, combined with the dif- 
ference in pressure probed by these wavelengths, indicates the presence of a temperature 
inversion. We find that a temperature inversion could occur somewhere between ~ 10 _1 - 
10~ 4 bar, where temperatures increase by roughly 500-700 K. This temperature-inverted 
region resembles those of planetary stratospheres; it occurs at a similar pressure level and 
causes a convectively stable temperature profile. The model spectrum from three exemplar 
temperature profiles (inverted) are shown in Fig. 3 (together with a temperature profile with 
no inversion), and the model residuals (Fig. 4) show the improved fit with a temperature 
inversion. The presence of the stratosphere coupled with a relatively steep stratospheric 
temperature gradient implies the existence of a local absorber that heats the atmosphere 
and maintains the local temperature struc ture. The temperature profile we derive con firms 
the previous identification of a T inversion teurrows. et al. 2007 ; Knutson et al. I2OO8I ). but 
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the profile reported here differs in terms of the detailed shape. 

Our estimates for the temperature profile and molecular mixing ratios are based on 
modeling the combined near-infrared and mid-infrared data sets. How ever, the recent finding 



of variability in the mid- infrared emission spectrum of HD 189733b (IGrillmair et al. 1120081 ) 
suggests caution is needed when deriving constraints from the noncontemporaneous mid-IR 
data. Although the presence of H 2 and CH 4 can explain most of the spectrum, we find that 
an additional absorber is required around 2.1 /im, and we attribute this to C0 2 . The evidence 
for a temperature inversion comes from the mid-infrared spectrum, and the mid-infrared data 
provide some additional constraints on the mixing ratios for H2O and CO2. We find that 
a range of temperature profiles and molecular mixing ratios are consistent with the data as 
illustrated by the models shown in Fig. 4 and the sample contribution functions shown in Fig. 
5, which correspond to three possible temperature profiles. Providing improved observational 
constraints on the location of the tropopause and the molecular mixing ratios will require 
additional observations, for which there are two possibilities. The first is to improve the 
spectral coverage and/or spectral resolution of the dayside spectrum, which could be done 
in the future with JWST. The second is to use a transmission spectrum (obtained during 
primary eclipse) to estimate the molecular abundances in the terminator regions; subject 
to some assumptions, this could be used to place constraints on the abundance of one or 
more molecules. HST transmission spectra for this planet have been obtained and are being 
analyzed by our team. However,the interpretation of the transmission spectrum is complex 
due to the presence of temporal variability (the subject of a forthcoming paper). 

There are two important caveats concerning our derived temperature profiles. First 
is the assumption of emission from a uniform disk rather than from a realistic irradiated 
hemisphere; as such, the spectrum and T-profile we derive is disk-averaged. Second, the 
spectral retrieval process used in our best fitting models implicitly acknowledges the possible 
role of dynamics in establishing atmospheric temperature by relaxing typical 1-line model 
constraints on the temperature gradient. Thus, portions of the atmosphere of HD 209458b 
may support departures from radiative equilibrium. Fully self-consistent modeling would 
require handling heat advection in the context of a global circulation model, which is beyond 
the scope of this paper. It is worth noting that the presence of a relatively strong dayside 
temperature inversion significantly complicates the spectral retrieval process for HD 209458b 
relative to the case of HD 189733b. 

Given that the dayside emission spectra for HD 189733b and HD 209458b have been 
observed with nearly identical instrument configurations, we can make a preliminary com- 
parison of these two planets (see Fig. 6). For the present, we restrict this discussion to 
the near-infrared to avoid the complications introduced by the presence of mid-infrared vari- 
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ability observed in HD 189733b (jGrillmair et al. II2008I ). In both planets, the near-infrared 



dayside emission spectrum probes the 5 < P < 10~ 2 bar portion of the atmosphere. Sim- 
ilarities - both planets show the presence of H 2 0, CO2, and 4p < for pressures near 1 
bar. Differences - the abundance of CH4 is significantly enhanced in HD 209458b relative to 
HD 189733b, and the temperature at 1 bar is higher in HD 209458b, as well. The observed 
differences in the near-infrared dayside spectra of these two hot-Jupiters are likely due to 
differences in temperature and composition; HD209458b is dominated by CH 4 absorption 
features, while HD 189733b is dominated by absorption from H2O and CO2. 



4. Conclusions 



In summary, we have presented the first near-infrared spectrum of light emitted by 
HD 209458b. Using an iterative forward model approach for spectral retrieval, we find that 
H 2 0, CH 4 , and CO2, together with a temperature inversion, are present in the dayside atmo- 
sphere of HD 209458b. There are a range of temperature profiles and molecular abundance 
solutions that are consistent with the data. Additional observational constraints on the at- 
mospheric temperature structure and composition will require either improved wavelength 
coverage/spectral resolution for the dayside spectrum or a transmission spectrum. We note 
that some of the temperature profiles consistent for these observations raise the question 
of whether the dayside atmosphere is in radiative equilibrium. Although advection of heat 
could support departures from radiative equilibrium, our present knowledge of most molec- 
ular opacities at high temperatures limits our ability to determine decisively whether the 
radiative equilibrium condition is met or not; thus there is an urgent need for further labo- 
ratory studies to obtain molecular databases for determining high temperature opacities of 
the most common molecules expected in hot- Jupiter atmospheres. 



We appreciate the Director's time award for these observations, and we thank Tommy 
Wiklind, Beth Padillo, and other members of the Space Telescope Science Institute staff for 
assistance in planning the observations. We also thank Jonathan Tennyson and Bob Barber 
for help with the water line list. G. Tinetti was supported by the Royal Society. A portion 
of the research described in this paper was carried out at the Jet Propulsion Laboratory, 
under a contact with the National Aeronautics and Space Administration. 
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Fig. 1. — This figure shows a set of normalized light curves plotted as a function of orbital 
phase. Gaps in the data are unavoidable because HD 209458 is not in the continuous 
viewing zone of Hubble. The top curve is the raw, broad-band, NICMOS data showing large 
internal scatter and a systematic offset affecting the first orbit. The second curve shows 
the calibrated, broad-band, NICMOS data (we exclude the first orbit); a shallow eclipse is 
visible. The third curve shows the calibrated 1.9 ^/m data. The larger eclipse depth at this 
wavelength (compared to the calibrated broad-band data) is readily apparent. The fourth 
curve shows the contemporaneous MOST visible photometry data. Both the individual and 
averaged data are shown. The MOST data are consistent with no changes in the system 
brightness at visible wavelengths. 
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Fig. 2. — A comparison of the observed spectrum (black points with ± l-a error bars) to a 
sequence of models showing the affect of H 2 0, CH 4 , and CO2. This comparison shows the 
portion of the near-infrared spectrum where each molecule makes a significant contribution. 
CH4 and H2O explain most of the spectral features, although additional absorption is needed 
around 2.0 jum (provided by CO2). Here, the observed spectrum has been smoothed with a 
5-point scrolling boxcar and thus adjacent data points are not statistically independent. 
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Fig. 3. — Top: Mid-infrared observations compared to synthetic spectra for three models 
with the tropopause at 0.001 bars (blue), 0.01 bars (purple), & 0.1 bars (green); these models 
illustrate the range of temperature/composition possibilities consistent with the data. We 
also show the best fitting model without a temperature inversion (red). Bottom: Model 
residuals compared to the mid-infrared observations. The model without a tropopause has 
the largest residuals. The contribution functions and temperature profiles associated with 
these models are shown in Fig. 5. The mid-infrared photometry are from Knutson et al. 
2008 and the mid-infrared spectroscopy are from Swain et al. 2008a. 
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Fig. 4. — The near-infrared and mid-infrared observations compared to synthetic spectra for 
three models (also shown in Fig. 3) that illustrate the range of temperature/composition 
possibilities consistent with the data. For each model case, the molecular abundance of 
CH 4 , H 2 0, & C0 2 and the location of the tropopause is given; the contribution functions 
for each of these models is shown in Fig. 5. The three dashed lines correspond to single- 
temperature thermal emission models with temperatures of 1400, 1800, and 2200 K; these 
serve to illustrate how the combination of molecular opacities and the temperature structure 
cause significant departures from a purely single-temperature thermal emission spectrum. 
Note that the mid-infrared data are not contemporaneous with the near-infrared data, and 
attempting to "connect" these data sets with a model spectrum is potentially problematic 
if significant variability is present. 
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Fig. 5. — The contribution functions for selected wavelengths and the associated temper- 
ature profile for each model shown in Figs. 3 & 4; starting from the top left and going 
clockwise, the models correspond to a tropopause at 0.001 bar, 0.01 bar, 0.1 bar, and no 
tropopause, respectively. The selected wavelengths correspond to bands of H 2 and CH 4 ; 
the contribution function is determined by opacity in each layer, which, in turn, can depend 
on the local temperature. 
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Fig. 6. — Top: The near-infrared dayside emission spectrum of HD 189733b and HD 209458b, 
showing the significant differences in the observed spectra. Middle: The brightness temper- 
ature spectrum of HD 189733b and HD 209458b. The spectra probe similar pressure scales 
in the dayside atmosphere. The difference in these spectra is due primarily to the presence 
of (1) significant CH 4 enhancement, and (2) increased temperature in HD 209458b relative 
to HD 189733b. Bottom: A preliminary comparison of HD 189733b and HD 209458b in 
terms of the temperature and molecular abundances. The abundance ranges shown here 
represent the range of solutions we currently identify as consistent with the data. Other 
temperature/abundance solutions may exist, and thus the results here should be interpreted 
as indicative. There is a suggestion of enhancement in the abundance of CH 4 and H 2 in 
HD 209458b relative to HD 189733b; this apparent enhancement needs to be confirmed via 
additional observations and detailed modeling. 



